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We analyze the implications of neutrino masses for the magnitude of neutrino magnetic moments. 
By considering electroweak radiative corrections to the neutrino mass, we derive model-independent 
naturalness upper bounds on neutrino magnetic moments, pi„, generated by physics above the 
electroweak scale. For Dirac neutrinos, the bound is several orders of magnitude more stringent than 
present experimental limits. However, for Majorana neutrinos the magnetic moment contribution 
to the mass is Yukawa suppressed. The bounds we derive for magnetic moments of Majorana 
neutrinos are weaker than present experimental limits if is generated by new physics at ~ 1 
TeV, and surpass current experimental sensitivity only for new physics scales > 10 - 100 TeV. The 
discovery of a neutrino magnetic moment near present limits would thus signify that neutrinos are 
Majorana particles. 



I. INTRODUCTION 

In the Standard Model (minimally extended to include 
non-zero neutrino mass) the neutrino magnetic moment 
is given by 
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An experimental observation of a magnetic moment 
larger than that given in Eq.Q would be an unequivo- 
cal indication of physics beyond the minimally extended 
Standard Model. Current laboratory limits are deter- 
mined via neutrino-electron scattering at low energies, 
with ii v < 1.5 x l(T 10 /i S and p, v < 0.7 x IO^Vb 
obtained from solar and reactor experiments, respec- 
tively. Slightly stronger bounds are obtained from astro- 
physics. Constraints on energy loss from astrophysical 
objects via the decay of plasmons into vv pairs restricts 
the neutrino magnetic moment to be /i„ < 3 x 10 -12 
Neutrino magnetic moments are reviewed in and 
recent work can be found in 0, [|| . 

It is possible to write down a naive relationship be- 
tween the size of ii v and m„. If a magnetic moment is 
generated by physics beyond the Standard Model (SM) 
at an energy scale A, as in Fig. [IJ,, we can generically 
express its value as 
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where e is the electric charge and G contains a combi- 
nation of coupling constants and loop factors. Removing 
the photon from the same diagram (Fig.^i) gives a con- 
tribution to the neutrino mass of order 



GA. 



(3) 



We thus have the relationship 
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which implies that it is difficult to simultaneously recon- 
cile a small neutrino mass and a large magnetic moment. 
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FIG. 1: a) Generic contribution to the neutrino magnetic mo- 
ment induced by physics beyond the standard model, b) Cor- 
responding contribution to the neutrino mass. The solid and 
wavy lines correspond to neutrinos and photons respectively, 
while the shaded circle denotes physics beyond the SM. 

However, it is well known that the naive restriction 
given in Eq.Q can be overcome via a careful choice 
for the new physics. For example, we may impose a 
symmetry to enforce = while allowing a non-zero 
value for \x v 0, 0, 0, 0] , or employ a spin suppres- 
sion mechanism to keep m„ small [lj]. Note however, 
that these symmetries are typically broken by Standard 
Model interactions. The original version of the well- 
known Voloshin mechanism involved imposing an 
SU(2) V symmetry, under which the left-handed neutrino 
and antincutrino (v and v c ) transform as a doublet. The 
Dirac mass term transforms as a triplet under this sym- 
metry and is thus forbidden, while the magnetic moment 
term is allowed as it transforms as a singlet. However, the 
SU(2) U symmetry is violated by SM gauge interactions. 
For Majorana neutrinos, the Voloshin mechanism may 
be implemented using flavor symmetries, such as those 
in |lll IT2I H^ . These flavor symmetries are not bro- 
ken by SM gauge interactions but are instead violated 
by SM Yukawa interactions, provided that the charged 
lcpton masses are generated via the standard mechanism 
through Yukawa couplings to the SM Higgs boson. 1 By 
calculating neutrino magnetic moment contributions to 



1 If the charged lepton masses are generated via a non-standard 
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m v generated by SM radiative corrections, we may thus 
obtain general, "naturalness" upper limits on the size of 
neutrino magnetic moments. 

In the case of Dirac neutrinos, a magnetic moment 
term will generically induce a radiative correction to the 
neutrino mass of order 15] 
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If we take A ~ 1 TeV and m„ < 0.3 eV, we obtain the 
limit \x v < 10~ 15 /x_b, which is several orders of magnitude 
more stringent than current experimental constraints. 

The case of Majorana neutrinos is more subtle, due to 
the relative flavor symmetries of m v and \i v respectively. 
The transition magnetic moment [A*i/]~g is antisymmetric 
in the flavor indices {a, j3}, while the mass terms [m v ] tt ^ 
are symmetric. These different flavor symmetries play 
an important role in our limits, and are the origin of the 
difference between the magnetic moment constraints for 
Dirac and Majorana neutrinos • 

It has been shown in that the constraints on Ma- 
jorana neutrinos are significantly weaker than those ob- 
tained in for Dirac neutrinos, as the different flavor 
symmetries of m„ and \x v lead to a mass term which is 
suppressed by charged lepton masses. This conclusion 
was reached by considering one-loop mixing of the mag- 
netic moment and mass operators generated by Standard 
Model interactions. The authors of Ref. found that 
if a magnetic moment arises through a coupling of the 
neutrinos to the neutral component of the SU(2)l gauge 
boson, the constraints obtained in [Tq for /j, Te and /j, TjJ- 
are comparable to present experiment limits, while the 
constraint on [i e{i is significantly weaker. Furthermore, 
constraints on magnetic moments generated only via a 
coupling to the U(l)y (hypercharge) gauge boson were 
found to be additionally suppressed. Thus, the analysis 
of |l6| lead to a model independent bound that is less 
stringent than present experimental limits. 

We shall show that two-loop matching of mass and 
magnetic moment operators implies stronger constraints 
than those obtained in if the scale of the new physics 
A > 10 TeV. Moreover, these constraints apply to a mag- 
netic moment generated by either the hypercharge or 
SU(2)l gauge boson. In arriving at these conclusions, 
we construct the most general set of operators that con- 
tribute at lowest order to the mass and magnetic mo- 
ments of Majorana neutrinos, and derive model indepen- 
dent constraints which link the two. We thus obtain a 
completely model independent naturalness bound that - 



mechanism, SM Yukawa interactions do not necessarily violate 
flavor symmetries. However, such flavor symmetries must always 
be broken via some mechanism in order to obtain non-degenerate 
masses for the charged leptons. 



for A > 100 TeV - is stronger than present experimen- 
tal limits (even for the weakest constrained element /x e/J .) 
Our key findings are summarized in Table [I] 

Our result implies that an experimental discovery of a 
magnetic moment near the present limits would signify 
(i) neutrinos are Majorana fermions and (ii) new lepton 
number violating physics responsible for the generation 
of \i v arises at a scale A which is well below the see-saw 
scale. Implications for neutrinoless double beta decay 
in theories with low scale lepton number violation have 
been discussed in [l7| . 



II. FRAMEWORK 



Following Refs. [15J, [16j we assume that the magnetic 
moment is generated by physics beyond the SM at an 
energy scale A above the electroweak scale. 2 In order to 
be completely model independent, the new physics will 
be left unspecified and we shall work exclusively with 
dimension D > 4 operators involving only SM fields, ob- 
tained by integrating out the physics above the scale A. 
We thus consider an effective theory that is valid below 
the scale A, respects the SU(2)l x U(1)y symmetry of 
the SM, and contains only SM fields charged under these 
gauge groups. We assume that all the usual SM interac- 
tions are present, including the Yukawa couplings of the 
charged leptons to the SM Higgs. We shall also work with 
the flavor states v a (i.e. the basis where the charged lep- 
ton masses are diagonal) and discuss the mass eigenstate 
basis in Section [V] Note that in either basis, Majorana 
neutrinos cannot have diagonal magnetic moments, but 
are permitted non-zero transition moments. 

The lowest order contribution to the neutrino (Ma- 
jorana) mass arises from the usual five dimensional op- 
erator containing Higgs and left-handed lepton doublet 
fields, L and H, respectively: 
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where e = — iT2, L c = L T C, C denotes charge conjuga- 
tion, and a, (3 are flavor indices. The neutrino magnetic 
moment operator must be generated by gauge invariant 
operators involving the SU(2)l and U(l)y gauge fields, 
W£ and -B M respectively. The lowest order contribution 
arises at dimension seven, so we consider the following 
operators, 

P B ] a0 = g' ' (p a eH) {H T eLp)B^, (7) 
PwU = 9 (^ eH ) a"" {H T er a L f3 ) W^, (8) 



,,. - dpBv - d v B„ and W°„ = d M IV„ a - d v W£ 



where B : 

ge a bcW*WZ are the U(l)y and SU(2) L field strength ten- 
sors, respectively, and g' and g are the corresponding 



2 We do not consider models with new light particles, such as Il8l . 
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couplings. We also define a 7D mass operator as 

[Ol?] ap = (Ts a eH)(H T eLp)(^H). (9) 

We note that the three 7D operators given by Eas. (|7l9|) 

do not form a closed set under renormalization. The full 
basis of 7D operators may be found in the appendix. 

Operators 0\f and 0\f are flavor symmetric, while 
Os is antisymmetric. The operator Ow is the most gen- 
eral 7D operator involving W® u (see the appendix for 
details). However, as it is neither flavor symmetric nor 
antisymmetric it is useful to express it in terms of opera- 
tors with explicit flavor symmetry, O^, which we define 
as 



[°w] aP = \{[Ow] a[j ±[O w ], a } 
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Our effective Lagrangian is therefore 
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■ ■ " denote other terms that are not rele- 
vant to the present analysis. 

Due to the differences in the flavor structure of O w 
and O w , the phenomenological manifestations of C w 
and C w are quite distinct. After spontaneous symme- 
try breaking, the flavor antisymmetric operators Ob and 
O w contribute to the magnetic moment interaction 

\ ( 12 ) 

where F^ v is the electromagnetic field strength tensor, 
W] a fi 2m e v 2 
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and the Higgs vacuum expectation value is (H T ) — 
(0,v/^/2). The flavor symmetric operator O w does not 
contribute to this interaction at tree-level. Similarly, the 
operators 0\f and Oj^ generate contributions to the 
Majorana neutrino mass terms, | [ m A a (3 vC uVpi given by 
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iii) 2-loop, 7D 

iv) 2-loop, 5D 
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TABLE I: Summary of constraints on the magnitude of the 
magnetic moment of Majorana neutrinos. The upper two 
lines correspond to a magnetic moment generated by the Ow 
operator, while the lower two lines correspond to the Og op- 
erator. 



III. SU(2) GAUGE BOSON 

We first consider contributions from O w to the mass 
operators. Contributions to 0\f occur through match- 
ing of the effective theory onto the full theory at the 
scale A, illustrated at one-loop order with the diagrams 
of Fig. |3 If the operator O^ is inserted at the vertex, 
these two diagrams sum to zero. This is no surprise, as 
these two diagrams contain no non-trivial flavor depen- 
dence. We therefore cannot obtain a contribution to the 
flavor symmetric mass operator by inserting the flavor 
antisymmetric operator O w . 





FIG. 2: Contribution of O w to the 5D neutrino mass op- 
erator. The shaded box indicates the O w operator while 
solid, dashed, and wavy lines correspond to leptons, Higgs, 
and gauge bosons, respectively. 



However, one-loop matching does yield a contribution 
to 0|P associated with the flavor symmetric operator 



O w of order 
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In sections lIHl and llVl below. we calculate radiative cor- 
rections to the neutrino mass operators (Ofp and 7 A f) 
generated by the magnetic moment operators O^ and 
Ob- This allows us to determine constraints on the size 
of the magnetic moment in terms of the neutrino mass, 
using Eas. (|13fl and Q14JI. Our results are summarized in 
Table [fl below, where the quantity R a p is defined as 
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with m a being the masses of charged lepton masses. Nu- 



merically, one has R Te ~ R T n ~ 1 and i? p 



283. 
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where we have estimated the one-loop matching contribu- 
tions using Fig. Eland naive dimensional analysis (NDA). 
We note that this contribution arises from loop momenta 
of order A, and that the precise numerical coefficient can- 
not be obtained without knowing the structure of the 
theory above this scale. We therefore follow the argu- 
ments of Rcf. 0| and employ NDA to estimate these 
high momentum contributions. 

We see that the one-loop contribution to the 5D mass 
term provides a strong constraint on C w but no con- 
straint on the parameter C w . In general, C^> are un- 
related parameters in the theory. If the new physics 
were to have no specific flavor symmetry/ antisymmetry 
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it might be natural for C^j> to be of similar magnitude. 
Alternatively, given the strong constraint on aris- 
ing from Eq. i|16|) . a sizable magnetic moment requires 
\C^y\ ^> \Cyy\. We identify two scenarios, to be discussed 
below: 



. \a 



in 



Arbitrary \C V 



A. Special case: C v 



r 4 



We have seen that the flavor antisymmetric operator 
0^ does not contribute to the 5D neutrino mass term at 
1-loop order; thus a direct constraint on the magnetic mo- 
ment is not obtained from the diagrams in Fig. [2] How- 
ever, suppose we had a theory in which the coefficients 
of O^y and were of similar magnitude, ~ C^. 
Then, using Eas.lfT5)l. (jT4"|) and (JTBJ) we have 
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and thus obtain a stringent \i v bound similar to that for 
Dirac neutrinos. 

We emphasize that Eq. (|17|l is not a model-independent 
constraint, as in general and O^r are unrelated. 
While it might seem natural for the the new physics to 
generate coefficients of similar size for both operators, we 
could obtain finite and vanishing (at tree-level) 
by imposing an appropriate flavor symmetry. 



where m a are the charged lepton masses, and the exact 
coefficient has been computed using dimensional regu- 
larization and renormalized with modified minimal sub- 
traction. We note that Eq. I|18ll gives the contribu- 
tion to the neutrino mass from all scales between A and 
the scale of electroweak symmetry breaking to leading 
log order. Moreover, while all contributions of the type 
(g 2 In A/Mw) n could be resummed using the renormal- 
ization group, a study of analogous operator mixing for 
Dirac neutrinos suggests that it is sufficient to retain 
only the leading log contribution. Using this result - as 
well as Eqs. l(T3)) and (fHEjl to relate and Cj^ to /z^ 
and m„ respectively - leads to bound (i) in Table [I] 

Note that this provides a weaker constraint than that 
in Ea. l|17[) . as it is suppressed by the charged lepton 
masses, and also only logarithmically dependent on the 
scale of new physics A. 




FIG. 3: Contribution of O w to the 7D neutrino mass opera- 
tor. 



2. 5D mass term — 0\ 



The neutrino magnetic moment operator O^y will also 
contribute to the 5D mass operator via two-loop match- 
ing of the effective theory onto the full theory at fj, ~ A. 



General case: Arbitrary C v 



We now consider the more general case where and 
C^, are unrelated, and directly derive constraints on the 
the coefficient of the flavor antisymmetric operator, C^. 




1. 7D mass term — Ow 

As the operator is flavor antisymmetric, it must 
be multiplied by another flavor antisymmetric contribu- 
tion in order to produce a flavor symmetric mass term. 
This can be accomplished through insertion of Yukawa 
couplings in the diagram shown in Fig. [3] 16). This di- 
agram provides a logarithmically divergent contribution 
to the 7D mass term, given by 



[C™(M w )] a 
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FIG. 4: Representative contribution of O w to the 5D neutrino 
mass operator. 

An illustrative contribution is shown in Fig.^J As with 
the diagrams in Fig. we require two Yukawa insertions 
in order to obtain a flavor symmetric result. This dia- 
gram contributes to the 5D mass operator, and we again 
provide an NDA estimate: 



af3 ' 



(19) 



(16tt 2 ) 2 v 2 

Again, using Eas. (|13fl and 114(1. this leads to bound (ii) 
m Table D In doing so, we have neglected the running 
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of the operator coefficients from the scale A to My since 
the effects are higher order in the gauge couplings and 
have a negligible numerical impact on our analysis. 

Compared to f-loop (7D) case of Eq. (|f 8fl . the 2-loop 
(5D) matching leads to a mass contribution that is sup- 
pressed by a factor of 1/16-7T 2 arising from the additional 
loop, but enhanced by a factor of A 2 jv 2 arising from the 
lower operator dimension. Thus, as we increase the new 
physics scale, A, this two-loop constraint rapidly becomes 
more restrictive and nominally provides a stronger con- 
straint than the f-loop result once A > Anv ~ 4 TeV. 
Inclusion of the logarithmic A-dependence of one-loop 
mixing implies that the "crossover " scale between the 
two effects is closer to ~ 10 TeV. 



IV. HYPERCHARGE GAUGE BOSON 

Unlike the case of the SU(2)l gauge boson, where a 
flavor symmetric operator O w exists, the operator Ob 
is purely flavor antisymmetric. Therefore, it cannot con- 
tribute to the 0\f mass term at one loop. As was noticed 
in (lfjj. the one- loop contribution of Ob to the O r A ® mass 
term also vanishes. 



3. 7D mass term — Ob 

If we insert Ob in the diagram in Fig. [3] the contribu- 
tion vanishes, due to the SU(2) structure of the graph. 
Therefore, to obtain a non-zero contribution to 0\f from 
Ob we require the presence of some non-trivial SU{2) 
structure. This can arise, for instance, from a virtual W 
boson loop as in Fig. |S] 16]. 




FIG. 5: Representative contribution of Ob to the 7D neutrino 
mass operator at two loop order. 

This mechanism gives the leading contribution of the 
operator Ob to the 7D mass term. The Ob and Ow 
contributions to the 7D mass term are thus related by 
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where 9w is the weak mixing angle and where the fac- 
tor on the RHS is due to the additional SU(2)l boson 
loop. This additional loop suppression for the Ob contri- 
bution results in a significantly weaker neutrino magnetic 



moment constraint than that obtained above O w . The 
corresponding limit is shown as bound (hi) in Table Q] 



4- 5D mass term — Ob 

However, the leading contribution of Ob to the 5D 
mass term arises from the same 2-loop matching consid- 
erations (Fig. 2J) that we discussed in connection with 
the O w operator. Therefore, the contribution to the 5D 
mass term is the same as that for Ow , except for a factor 
of (g'/g) 2 = tan 2 6w- We thus obtain 
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corresponding to bound (iv) in Table. [U Importantly, 
this is the strongest constraint on the Ob contribution 
to the neutrino magnetic moment for any value of A. 



FLAVOR STRUCTURE 



All the magnetic moment bounds presented in Table [I] 
contain the factor R a p = m 2 /{m 2 a — rn 2 ^). This factor 
arises from the Yukawa insertions that are necessary to 
obtain flavor symmetric mass terms from the flavor an- 
tisymmetric magnetic moment operators. The strongest 
constraints are those for the ji re and fi Tfi transition mo- 
ments, as R Te ~ i? TM ~ 1. However i?^ e = 283, so the 
constraint on the e—fi transition moment is much weaker. 

In the foregoing analysis, it was convenient to work in 
the flavor basis (where the charged lepton masses are di- 
agonal). However, given /z Q| g, we can use neutrino mixing 
matrix to determine fj,^, where i,j denote neutrino mass 
eigenstates. Rotating the neutrino magnetic moments 
from the flavor basis to the mass basis gives 



a/3 



(22) 



Since most elements of Uij are large, it would be natural 
to expect that all elements of would be of similar size. 
In particular, all elements of will receive a contribu- 
tion from (the most weakly constrained element in 
the flavor basis.) Therefore the constraint will trans- 
late into similar limits for all elements of . 



VI. DISCUSSION AND CONCLUSIONS 

We have discussed radiative corrections to the neutrino 
mass arising from a neutrino magnetic moment coupling. 
Expressing the magnetic moment in terms of effective 
operators in a model independent fashion required con- 
structing operators containing the SU(2)l and hyper- 
charge gauge bosons, Ow and Ob respectively, rather 
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than working directly with the electromagnetic gauge bo- 
son. We then calculated //„ naturalness bounds arising 
from the leading order contributions to both the 5D and 
7D Majorana mass terms. These bounds are summarized 
in Table IJ 

At the TeV scale, the strongest bound comes from a 1- 
loop contribution to the 7D mass term as previously dis- 
cussed in [l6( [limit (i) in Table Q]. However, this bound 
applies only to a /ij; contribution to the magnetic mo- 
ment, and not to a contribution. Even for the ^ 
contribution, the bound on the e — ji transition moment 
is weaker than current experimental limits. 

However, we have shown that 2-loop contributions to 
the 5D mass term always provide the strongest bound on 
jti^ [limit (iv) in Table d and also provides the strongest 
bound on provided A > 10 TeV [limit (ii) in Tabled. 
We have estimated this contribution - associated with 
matching the effective theory onto the (unspecified) full 
theory at high scales - using NDA. The resulting expres- 
sions for and ^ differ only by a factor of tan 2 6w 
Taking the [i„ limit (the weaker of the two) we have 

s a m-9 fKL/sN ( 1 TeV\ 2 m 2 T 

(23) 

For any scale A, Ea . (|23|) is the most general naturalness 
bound on the size of the Majorana neutrino magnetic 
moment. It can only be evaded in the presence of fine 
tuning of the couplings or model-dependent suppression 
of the matching conditions at the scale A. 

Turning now to the current experimental situation, the 
best laboratory limit obtained from scattering of low en- 
ergy reactor neutrinos is "fjt, e " < 0.7 x 10~ 10 /U_b 0. Note 
that this limit applies to both and ^ Me , as the flavor 
of the scattered neutrino is not detected in the experi- 
ment. Taking the neutrino mass to be m v i^0.3 eV (as 
implied by cosmological observations, e.g.[2£|), Eq.J23) 
gives 

Mr^Mre < 1(T 9 [A(TeV)] ~ 2 

(i,e<3x 10~ 7 A(TcV)p 2 . (24) 

We thus conclude that if fi^e is dominant over the 
other flavor elements, an experimental discovery near the 
present limits (e.g., at fj, ~ lO^ 11 ^) would imply that 
A < 100 TeV. On the other hand, any model that leads 
to all elements of having similar size would only be 
consistent with the m v naturalness bounds if A < 10 
TeV. 

Our conclusions can be summarized according to the 
scale of the new physics, A: 

I) A < 10 TeV 

• No conflict with experimental limits. 

• Both 0\y and Ob contributions to [i v are 
possible, though Ow contributions are more 
tightly constrained. 



II) A > 10 TeV 

• jU T /ij Mre bounds stronger than experimental limits 

• jifie bound weaker than experimental limits 

• Same limit irrespective of whether [l v generated 
by Ow and Ob 

III) A > 100 TeV 

• The n a /3 bound becomes stronger than current 
experimental constraints, for all a, [3. 

The limit on the the magnetic moment of a Dirac neu- 
trino is considerably more stringent than for Majorana 
neutrino. This is due to the different flavor symmetries 
involved. In the Dirac case, no insertion of Yukawa cou- 
plings is needed to convert a flavor antisymmetric oper- 
ator into a flavor symmetric operator, and the stringent 
limit \x < 10 _15 yLts holds (in the absence of strong can- 
cellations) . A significant implication is that if a magnetic 
moment /i > 10~ 15 /xs were measured, it would indicate 
that neutrinos are Majorana fermions, rather than Dirac. 
Moreover, the scale of lepton number violation would be 
well below the conventional see-saw scale. 
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APPENDIX A: FULL 7 DIMENSIONAL 
OPERATOR BASIS FOR MAJORANA 
NEUTRINOS 

Here, we analyze the full set of gauge invariant oper- 
ators of dimensions five through seven containing left- 
handed lepton doublet fields, one Higgs doublet, and 
SU(2)l and U(1)y gauge bosons. It is convenient to 
classify the operators according to their dimension and 
field content. 

Mass operator (two fermion and four Higgs fields): 

7 A f = (L°eH) (H T eL) (H*H) (Al) 

Magnetic moment operators (two fermions, two Higgs 
and one gauge boson field): 

B = (L c eH)a^(H T eL)B^ 

O a w = (L-eH)a^(H T er a L)wr 

O b w = (L-eT a H)<j^(H T eL)wr 

O c w = (L-ea^L)(HeT a H)Wr (A2) 
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However, it can be shown that the operators involving 
W£ v are related via 

20^ -O^ = -20^ + 0^ = O a w -O h w . (A3) 

Furthermore, we have 

[°w] aP = [OfrW (A4) 

Therefore, there is only one independent operator involv- 
ing W£ v , which we may choose to be Qm = Ow- We 
note that the operator defined in Ref. [l(|, 

Davidson = ^er a a^L) (HeT b H) W^e abd , (A5) 

corresponds to the flavor antisymmetric component of 
O w . Explicitly, O° avidson = O w -O b w = 20 w . 

Two fermion fields, two Higgs fields, and two derivatives: 

Of = {L c eH) (DuJE^eD^L) 
0\ = (D^eD^H) (H T eL) 
Of = (D^eH) (H T eD»L) 
Of = (L c eD^H) (D^H T eL) (A6) 

It is also possible to establish relationships between some 
of the Of operators. Transposition yields: 

(Off = - (D (1 L T e T D^H) (H T e T C T L) 

= Ol (A7) 

where the minus sign arises from the interchange of two 
fermion fields, and we use e T = — e, C T = —C. Since 
the transpose of a number is the same number again, we 
conclude that there is only one independent operator out 
of these two, which we may choose to be Of = 0\. 

The same trick applied to the remaining two operators 
does not lead to any new relations. However, another 
relation can be obtained via integration by parts. We 
consider 

= D M (D^L c eHH T eL) = {11} + {12} + {13} + {14}, 

where the numbers in square brackets denote the field on 
which the first and the second derivative act, respectively. 
The squared derivative acting on Higgs is zero by the 



equation of motion, and that acting on fermion fields 
leads to operators of the form Ow,b- Repeating this trick 
for all possible combinations, we obtain a system of 4 
equations: 

{11} + {12} + {13} + {14} = (A9) 
{12} + {23} + {24} = 
{13} + {23} + {34} = 

{14} + {24} + {34} + {44} = 

i({ll} + {44}) + {12} + {13} 

+{14} + {23} + {24} + {34} = 0, 
where the last equation is the sum of the first four. Sub- 
tracting from this 5th equation the 1st and the 4th, we 
obtain: 

{14} = {23}-i({ll} + {44}), (A10) 

which implies that the operator Of can be expressed in 
terms of the operators O3 and Ow,b- We are thus left 
with just two independent operators with two derivatives: 

0\ = (L c eH) (D^eD^L) 

Ol = (L c eD^H) (D^H T eL) . (All) 

Two fermion fields, three Higgs fields, and one derivative: 

Ol = e c 7 M (H T er a L) (D^H T er a H) 

0\ = (D»H T eT a L) (H T 6T a H) (A12) 

After fierzing, both operators reduce to just one inde- 
pendent operator: 

4 = e c 7 M (H T eL) (D»H T eH) (A13) 

Four fermion fields and one Higgs: 

Ol = (L c eL) (Lje c ) ,,,11, 

Oj = {L^a, v L){Lja^e c )e l0 H 3 (A14) 

Thus {0\f, O b , Ow, 0|, Of , O4, Of , O^} form the 
complete set of 7D operators (where flavor labels have 
been suppressed). 



[1] W. J. Marciano and A. I. Sanda, Phys. Lett. B 67, 303 
(1977); B. W. Lee and R. E. Shrock, Phys. Rev. D 16, 
1444 (1977); K. Fujikawa and R. Shrock, Phys. Rev. Lett. 
45, 963 (1980). 

[2] J. F. Beacom and P. Vogel, Phys. Rev. Lett. 83, 5222 
(1999); D. W. Liu et ai, Phys. Rev. Lett. 93, 021802 
(2004). 

[3] H. T. Wong [TEXONO Collaboration], hep-ex/0605006 
B. Xin et al. [TEXONO Collaboration], Phys. Rev. D 72, 



012006 (2005); Z. Daraktchieva et al. [MUNU Collabora- 
tion], Phys. Lett. B 615, 153 (2005). 

[4] G.G. Raffelt, Phys. Rep. 320, 319 (1999). 

[5] M. Fukugita and T. Yanagida, Physics of neutrinos 
and applications to astrophysics, Chapter 10, Springer, 
Berlin, (2003), and references therein. 

[6] B. Kayser, F. Gibrat-Debu and F. Perrier, World Sci. 
Lect. Notes Phys. 25, 1 (1989). 

[7] H. T. Wong and H. B. Li, Mod. Phys. Lett. A 20, 1103 



8 



(2005). 

[8] G. C. McLaughlin and J. N. Ng, Phys. Lett. B 470, 157 
(1999); R. N. Mohapatra, S. P. Ng and H. b. Yu, Phys. 
Rev. D 70, 057301 (2 004). 
[9] A. B. Balantekin, |hep-ph/060TTT3| A. de Gou- 
vea and J. Jenkins, |hep -ph /0603036 A. Friedland, 
hep-ph/0505165 

[10] M. B. Voloshin, Sov. J. Nucl. Phys. 48, 512 (1988). For 
a specific implementation, see R. Barbieri and R. N. Mo- 
hapatra, Phys. Lett. B 218, 225 (1989). 

[11] H. Georgi and L. Randall, Phys. Lett. B 244, 196 (1990). 

[12] W. Grimus and H. Neufeld, Nucl. Phys. B 351, 115 
(1991). 

[13] K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 64, 
1705 (1990). 



[14] S. M. Barr, E. M. Freire and A. Zee, Phys. Rev. Lett. 

65, 2626 (1990). 
[15] N. F. Bell, V. Cirigliano, M. J. Ramsey-Musolf, P. Vogel 

and M. B. Wise, Phys. Rev. Lett. 95, 151802 (2005). 
[16] S. Davidson, M. Gorbahn and A. Santamaria, Phys. Lett. 

B 626, 151 (2005). 
[17] V. Cirigliano, A. Kurylov, M. J. Ramsey-Musolf and 

P. Vogel, Phys. Rev. Lett. 93, 231802 (2004). 
[18] K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 63, 

228 (1989). 

[19] C. P. Burgess and D. London, Phys. Rev. D 48, 4337 
(1993). 

[20] D. N. Spergel et al, astro-ph/0603449 



